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Merger induced scatter and bias in the cluster mass 
Sunyaev-Zeldovich effect scaling relation 
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ABSTRACT 

We examine sources of scatter in scaling relations between galaxy cluster mass and thermal 
Sunyaev-Zeldovich (SZ) effect using cluster samples extracted from cosmological hydrody- 
namical simulations. Overall, the scatter of the mass-SZ scaling relation is well correlated 
with the scatter in the mass-halo concentration relation with more concentrated halos having 
stronger integrated SZ signals at fixed mass. Additional sources of intrinsic scatter are projec- 
tion effects from correlated structures, which cause the distribution of scatter to deviate from 
log-normality and skew it towards higher inferred masses, and the dynamical state of clusters. 
We study the evolution of merging clusters based on simulations of 39 clusters and their cos- 
mological environment with high time resolution. This sample enables us to study for the first 
time the detailed evolution of merging clusters around the scaling relation for a cosmologi- 
cally representative distribution of merger parameters. Major mergers cause an asymmetric 
scatter such that the inferred mass of merging systems is biased low. We find mergers to be 
the dominant source of bias towards low inferred masses: over 50% of outliers on this side 
of the scaling relation underwent a major merger within the last Gigayear. As the fraction 
of dynamically disturbed clusters increases with redshift, our analysis indicates that mergers 
cause a redshift-dependent bias in scaling relations. Furthermore, we find the SZ morphology 
of massive clusters to be well correlated with the clusters' dynamical state, suggesting that 
morphology may be used to constrain merger fractions and identify merger-induced outliers 
of the scaling relation. 



Key words: cosmology: large-scale structure of Universe, galaxies: clusters: general, meth- 
ods: N-body simulations 



1 INTRODUCTION 

Clusters of galaxies are the most massive gravitationally bound 
objects in the universe, which makes them an important tool for 
cosmology: among other tests, their abundance provides informa- 
tion on the gravitational growth of structures and is regulated by 
the initial density field, gravity, and the expansion history of the 
universe, which critically depend on the underlying cosmology. 
Thus number counts of clusters, for which masses and redshifts 
are known, can b e used to constrain cosmological parameters (see 
I Allen et alJl201lL for a recent review). 

To relate observed number counts to theoretical predictions of 
the cluster mass function, these experiments need to infer cluster 
masses from observables. The thermal Sunyaev Zeldovich (SZ) 
efl^ect, the signature of inverse Compton scattering of cosmic 
microwave background photons with hot cluster electron, is 
thought to provide an excellent mass proxy as the SZ signal is 
proportional to the total thermal energy of a cluster and is thus less 



afi^ected by physical processes in the cluster core which can largely 
affect the X-ray luminosity. This is confirmed by simulations (e, 
iNagai 2006; Shaw et al. 2008 : Battaglia et al. 2010 ; Sehgal et 
|2010) finding the scatter in the mass - SZ scaling relation to be of 
order 5 - 10%. Furthermore, the SZ effect is not subject to surface 
brightness dimming and has a very weak redshift dependence, 
making it an ideal probe to study galaxy clusters at high redshift. 



Currently several large surveys are sta rting to detect hundreds 
of galaxy clusters through their SZ signal ([V anderlind e et al J20ld ; 
[Marriage et alj2010l ; |Planck Collaboration et al.201 la[) and derive 



cosm o logical constraint ba s ed on these samp les i Andersson et al.l 
l20ld ; ISehgal et al.l |2010| ; IWilliamson etaU l201lh . To exploit 
the statistical power of these upcoming cluster samples, the 
mapping between SZ signal and cluster mass needs to be well 
understood. Observations find normalization and slope of the 
scaling relations bet ween SZ signal and lensing derived masses 
JMarrone et alj|201 ih . or between SZ signal and X-ray properties 
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jPlanck Collaboration et alj 1201 Ibl ld) to be consistent with self- 
similar scaling and predictions from simulations. 

Due to the steep slope of the cluster mass function, com- 
petitive cosmological constraints from these experiments require 
information about the distribution a nd redshift evolution of 
scatter in the mass scaling relation (e.g. iMaiumdar & Mohnl2004l : 
iLima & Hij|2005l : IShaw et al.ll201(]h . As the true cluster mass and 
other physical cluster properties which may bias the mass proxy 
are unobservable, and as the noise and biases in the different mass 
estimators may be correlated, characterizing the intrinsic scatter in 
any of these scaling relation is difficult to obtain from observations. 
Hence the sources and distribution of scatter in different mass 
estimators are mainly studied through simulations and mock 
obser v ations (e.g.lRasiaet alj2006:lNagai et al.l2007q:IShaw et al 



20111) 



2004 iBecker & Kravtsovl I2OI0I : lYang et al.1 I2OI0I ; fpabian et al 



In this work we focus on the effect of merging events on the 
SZ signal of a galaxy cluster. As clusters form through merging 
of smaller objects these are frequent and disruptive events, which 
may alter the physical state of the involved clusters significantly. 
Hence merging clusters may deviate from the scaling relations 
observed in relaxed clusters and, as the fraction of morphologi- 
cally disturbed clusters increases with redshift, cause a redshift 
dependent scatter or bias in the mass scaling relat i on. Simula- 
tions of bin a ry cluster mergers JRandall et al.l I2OO2I : IPoole et al.l 
I2OO6I l2007l : IWik et alj l2008h find that the X-ray luminosities, 
temperatures, SZ central Compton parameters and integrated SZ 
fluxes increase rapidly during the first and second passage of 
the merging clusters. The clusters temporarily drift away from 
mass scaling relations and return to their initial scaling relation 
as the merging system virializes. These transi ent merger boosts 
found in binary mergers and some observations dSmith et alj2003r) 
can scatter the inferred masses towards higher values and thus 
bias the derived cosmology towards a higher nor malization of 
the power spectrum, a^ , and lower matter density JRandall et al.l 
200i ISmith et all |2003| ; IWik et al.l I2OO8I : lAngrick & Bartelmamj 
201 ih . On the other hand, mergers increase the non-thermal 



press ure support JRasia et al.ll2006l : lLau et al. 2009; Battaglia et al. 
I2OIOI) found in cluster outskirts, and due to partial virialization 
merging cluster s can appear cooler than r elaxed clusters of the 
same mass (e.g. lMathiesen & Evrardll200lh. For a cluster sample 
extracted from cosmological simulations, iKravtsov et alj <2006l) 
find the X-ray temperatures of morphologically disturbed clusters 
to be biased, while the X-ray derived SZ-equivalent Yx shows no 
significant correlation with cluster structur e. Comparing X-ray 
and SZ to weak lensing derived masses, lOkabe et al.l j2010f) 
and iMarrone et al.l ( 120111) found undisturbed clusters to have of 
order ~ 40% higher weak lensing masses than disturbed clusters at 
fixed T and Ksz, and ~ 20% higher weak lensing masses at fixed Kx- 

Our goal is to isolate how mergers in a cosmological context 
affect the SZ signal of clusters, and if merging cluster can be de- 
tected as outliers of scaling relations. This extends previous work, 
as our analysis includes both multiple mergers with realistic dis- 
tributions of orbits and mass ratios, and full SPH treatment of gas 
physics with radiative cooling, star formation and supernova feed- 
back. The simulations and the cluster sample are described in Sect. 
[2I We discussion the best-fit scaling relations and their scatter in 
Sect. [3] The effect of merging events of the clusters SZ signal is 
quantified and the evolution of merging clusters with respect to the 



scaling relations is discussed in Sect.|4l In Sect.[5]we investigate if 
the dynamical state of clusters can be inferred from the morphol- 
ogy of the SZ signal. We summarize our results and conclude in 
Sect.m 



2 SIMULATIONS 

This analysis is based on two samples of galaxy clusters extracted 
from cosmological hydrodynamics simulations. In this section we 
summarize the simulated physics and describe the derived quanti- 
ties used in our analysis. 



2.1 Cluster samples 

Sample A To study the time evolution of the cluster SZ signal 
we use a sample of 39 galaxy groups and clusters with virial 

masses abo ve 3 x lO' ^Mg/Zi from simulations presented in 

iDolag et al.l (200g, l2009f) . 25 of these clusters are more massive 
than 10'''Mo//i. These structures were identified as 10 different 
regio ns in a (479Mpc//t)^ dark-matter-only cosmological simula- 
tion jYoshida et al.l uOOlh . and re-simulated at higher resolutio n 
using the Zoomed Initial Conditions method ^Tormen et al.lll997). 
The re-simulations, described in detail in Dolag et al. ( 20060 . 
are carried out with GADGET-2 ( ISpringelll2005l) . and include a 
uniform, evolving UV-background and radiative cooling assuming 
an optically thin gas of primordial composition. Star formation is 
included u sing the two-phase m o del of the interstellar medium 
(ISM) by ISpringel & Hemguistl J2003h . In this sub-resolution 
model the ISM is described as cold clouds, providing a reservoir 
for star formation, embedded in the hot phase of the ISM. Star for- 
mation is self-regulated through energy injection from supemovae 
evaporating the cold phase. Additional feedback is incorporated in 
the form of galactic winds triggered b y supernovae that drive mass 
outflows dSpringel & Hemauistll2003l) . 

The simulation assumes a flat ACDM cosmology with 
(nm,nb,o"8,^) = (0.3,0.04,0.9,0.7). It has a mass resolution 
of m-DM = 1.1 X IO'^Mq/H and mgas = 1.7 x IO^Mq/H and the 
physical softening length is e = 5kpc//i over the redshift range 
of interest. Our analysis is based on 52 snapshots covering the 
redshift range z = 1 to z = and separated evenly in time with a 
spacing of 154 Myrs between snapshots. 



Sample B The second cluster sampl e is a volume-l i mited sample 
of 117 clusters at z = described in lBorgani et alj ([2004). These 
clusters are identified in a (192Mpc//!)^ cosmological SPH simu- 
lation carried out with GADGET-2 and using the same physics as 
described above. This simulation assumes a flat ACDM cosmology 
with (n„„f2b,o-8,/i) = (0.3,0.04,0.8,0.7). The mass resolution is 
'"DM = 4.6 X 10'Mo//i and nig^s = 6.9 x 10*Mo//i, the physical 
softening length at z = is e = 7.5kpc//i. 



2.2 Masses and merging histories 

Halos are identified using a friend-of-friends algorithm and the 
cluster center is defined by the particle in a halo with the minimum 
gravitational potential. Cluster radii R^ and masses Ma are defined 
through spherical regions around the cluster center within which 
the average density is A times the critical density of the universe 



Jo 



p{r)Anr'dr= — i^^Ap, 



Ma 



(1) 



Merger induced scatter in the M-Y relation 3 



2.5 
2 2.0 
I 1.5 
"° 1.0 

0.3 
I 0.2 

CD 

•*-^ 0.1 
0.0 



all 

M>10"M3jh 


- 




1 
\ 


\ 


- 






\ 

1 


\ 
\ 




,- 




J-- 


i 




_i 


- 


~i — I- 


_J 




^--|.— ., 1 









0.5 



0.6 



0.7 0.8 

a 



0.9 



1.0 



Figure 1. Fractions accretion rate per unit redsliift (top panel) and merger 
fraction as a function of scale factor The solid shows the complete sample 
A, the dash-dotted line a subsample of massive clusters. The dotted line 
indicates the overall mean accretion rate. Accretion rate (merger fraction) 
are averaged over 3 (5) neighboring simulation snapshots to reduce noise. 



We identify mergers by a mass jump criteria applied to the 
mass history of the main progenitor. Motivated by the findings 
that the average mass accretion history of halos is well de- 
scribed by exponent ial growth with redshift JWechsler et al .11 20021 : 
iMcBride et al.l l2009) and that the average merger rate per halo per 
unit redshift is nearly constant for a wide range of halo masses 
and redshifts iFakhouri & Mai 120081) , we select merging events 
based on a threshold in fractional mass accretion rate per unit 
redshift dM/dz/M > f,^. We choose ^^ such that halos accrete on 
average 30% of the mass accreted since its formation redshift Zf, 
defined as the redshift at which a halo reaches half its present day 
mass, during mergers. We checked that our results are insensitive 
to the exact choice of f^: We find similar trends for any merger 
definition ^m > {dM/dz/M)^,au!,ta that requires the accretion rate 
dM/dz/M during mergers to be larger than the mean accretion rate 
(cf. discussion of Fig. [6}. 



in their re-simulation region live in denser regions than an aver- 
age cluster of the same mass in a volume limited sample. In the 
following discussion we refer to the most massive objects in their 
respective re-simualtion region as primary clusters, and all others 
as secondary clusters. 

Simulations indicate a dependence of halo formation histories 
on environment with merger being more frequent in dense environ- 
ments and late-forming massive clusters living in dense r environ- 
ment s than earlier forming cl usters of the same ma ss JGao et al.l 
I2OO5I : IWechsler et all I2OO6I : Ipakhouri & m3 l2009h . Hence the 
merging histories of cluster sample A might not be representative of 
those of a volume limited sample. To assess the impact of our sam- 
ple selection on halo formation histories we compare the formation 
redshifts of primary and seconda ry clusters in sample A and halos 
in the Millenium run simulation l lSpringel et alj2005h in Fig.|2] 

The symbols show the present day masses and formation red- 
shift Zf for all clusters in sample A. Primary clusters are indicated 
by star symbols. The dashed and dotted lines are a fit to the mean 
form ation time and i ts Icr scatter for halos in the Millenium Run 
from iMcBride et alj ( |2009[) . We convert the fitting formula from 
friend-of-friends halo mass with linking length b = 0.2 to Af2oo as- 
suming a constant conversion factor M200 = 0.7MpoF- For the mass 
range of our sample this conversion underestimates M200LI and bi- 
ases the fit for Zf to more recent times. 

Due to the differences in matter density used in simulation A 
(CIm = 0.3) and in the Millenium Run (Q.^ = 0.25) the average 
clusters in simulation A forms earlier than a cluster of the same 
mass in the Millenium Run. Hence formation redshifts for primary 
clusters in sample A are broadly consistent with the formation his- 
tory of halos in the Millenium run. Figure|2]suggests that secondary 
clusters in sample A may form somewhat later than primary clus- 
ters of the same mass. However, the distribution of formation red- 
shifts at fixed mass is not expected to be symmetric but to have a 
long tail towards later formation times and the comparison is lim- 
ited by the small number objects. Overall, we expect the merging 
histories analyzed in this study to be similar to those found in a 
volume limited sample. 



Figure[T]confirms that this merger definition does not strongly 
depend on cluster mass or redshift. The top panel shows the mean 
accretion rate as a function of scale factor for all clusters (solid 
line) and massive clusters (M > IO'^'Mq/A, dash dotted line), and 
the overall mean accretion rate (dotted line). The lower panel shows 
the fraction of clusters that are merging as a function of scale factor. 
There is a peak of merging activity around a = 0.9, but the accretion 
rate and merger fraction show no clear trends with cluster mass or 
redshift. 



2.2.1 Comparison to the Millenium Run 

The 39 cluster and group-scale sized halos in sample A are ex- 
tracted from 10 re-simulation regions selected from a large simula- 
tion box. One of the re-simulated regions hosts a filamentary struc- 
ture with four massive clusters (M > IO'^Mq/Zi), and three of the 
re-simulation regions hosting other massive clusters contain several 
other smaller clusters. The re-simulation technique allows us to an- 
alyze the evolution of these regions of interest in their cosmological 
context at a higher resolution. As a result of the re-simulation strat- 
egy, the mass distribution of this sample does not follow the cluster 
mass function, and clusters which are not the most massive object 



2.3 SZ maps 

The amplitude of the thermal SZ effect along a line of sight is pro- 
portional to the Compton y parameter 



meC^ 



I din. 



(2) 



where We and T, are the electron density and temperature, ^b is the 
Boltzmann constant, ctj the Thomson cross section, m, the elec- 
tron rest mass, and c the speed of light. For each cluster we ana- 
lyze Compton y parameter maps obtained from three orthogonal 
lines of sight. For sample A the projection depth is 8 Mpc and 
maps are produced using the map making tool Smac JDolag et al.l 



' For equal mass particles, a FOF grou p vyith linkin g length b is bounded 
by a surface of density Sil^ pait/('2.7Tb^) IWhitel2002h . Assuming that halos 
follow NFW-profiles with concentration c = (4, 7, 10), the ration between 
M200 and MpoF with b = 0.2 in the Millenium run cosmology is given by 
(0.71,0.80,0.85). In practice however, the conversion between these mass 
definitions is complicated by deviations from the NFW-profile and spherical 
symmetry. 
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Figure 2. Distribution of formation redshifts ^f as afunction of cluster mass. 
Tlie symbols represent clusters from sample A, the most massive cluster in 
each re-simulation region is marked with a star symbol. Dashed (dotted) 
lines show a fit to the mean (Icr scatter) formation redshift as a function o f 
friends-of-fiiends mass found in the Millenium Run I McBride et alj|2009l) . 
converted to spherical overdensity mass assuming M200 = O.VMpoF (see 
text for details). Formation redshift is defined as the redshift at which a halo 
reaches half its present day mass. One cluster in sample A forms before 
z = 1 , indicated by the left an'ow. Open circles indicate the clusters shown 
as examples in subsequent plots , labels indicate the names of these clusters 
in Table 1 of lDolag et all J2009l) . 



l2005f) and the JobRunner web applicatior|j. For sample B we use 
projected map s whic h include all material with 6i?vii described in 
lAmeplio et al.l (20071). From these maps we measure integrated Ka 
parameters within different overdensity radii (i?2500,^500,^200,^vir) 



Ya = 









dV n.T. 



(3) 



where the integration volume is a cylinder of radius Rt^ and height 
8 Mpc (or 12 i?vii) for sample A (or B). This definition of the inte- 
grated Y parameter includes projection effects due to halo triaxiality 
and nearby structures within the projection cylinder, but does not 
account for projection effects from uncorrelated large scale struc- 
ture along the line of sight. 



3 MASS SCALING RELATIONS 

Self-similar clusters models predict the gas temperature to scale as 

T oc {ME{z))-'^ . (4) 

Hence the self-similar prediction for the relation between integrated 
Compton Y parameter and mass is 



Fa oc M,^,,,t,T oc f^.MfE'^Hz) 



(5) 



In this section we determine the best fit scaling relations for the sim- 
ulated clusters and discuss the scatter in these relations, focussing 
on the role of mergers. 



Table 1. Best fit Ma(Ja) scaling relation parameters (Eq. (6))and logarith- 
mic scatter <tm at fixed Y, defined analogously to Eq.[8] A*/B* denote sam- 
ples A/B restricted to clusters at z = with M > 2 X IO'^Mq/Zi. 



Sample 


A 


A{z = 0) 


a 


P 


O-M 


A 


200 


-0.348 ± 0.007 


0.639 ±0.010 


-0.57 ± 0.08 


0.063 


A* 


200 


-0.281 ±0.042 


0.588 ± 0.020 


- 


0.042 


B 


200 


-0.297 ± 0.006 


0.617 ±0.007 


- 


0.042 


B* 


200 


-0.261 ±0.014 


0.593 ±0.010 


- 


0.027 


A 


500 


-0.466 ± 0.001 


0.641 ± 0.007 


-0.74 ±0.10 


0.089 


A* 


500 


-0.406 ± 0.036 


0.607 ± 0.020 


- 


0.042 


B 


500 


-0.400 ± 0.004 


0.626 ± 0.005 


- 


0.037 


B* 


500 


-0.379 ±0.011 


0.604 ± 0.009 


- 


0.024 



Table 2. Best fit Yi^(Mi^) scaling relation parameters (Eq. (7)) and loga- 
rithmic scatter cry at fixed mass. A*/B* denote sample A/B restricted to 
clusters at z = with M > 2 X IO'^'Mq/Zj. 



Sample 


A 


B(z = 0) 


y 


S 


a-y 


A 


200 


0.547 ± 0.003 


1.560 ±0.014 


0.85 ±0.10 


0.103 


A* 


200 


0.489 ± 0.052 


1.648 ±0.056 


- 


0.070 


B 


200 


0.494 ± 0.005 


1.555 ±0.017 


- 


0.071 


B* 


200 


0.445 ± 0.030 


1.668 ±0.044 


- 


0.046 


A 


500 


0.714 ±0.003 


1.553 ±0.017 


1.03 ±0.14 


0.136 


A* 


500 


0.697 ± 0.038 


1.601 ±0.051 


- 


0.068 


B 


500 


0.641 ± 0.003 


1.556 ±0.014 


- 


0.059 


B* 


500 


0.624 ± 0.013 


1.637 ±0.027 


- 


0.037 



3.1 Best fit scaling relations 

We now determine the best fit Ma(Ka) scaling relation 

(kpc" 
and KaC^^a) scaling relation 



Ma(Ka) = 10" 



E^iz) lO'^Me/h 



rA(MA) = 10" 



Ma 



10'-*Mo//i 



E\z) kpc^ 



(6) 



(7) 



where the self similar predictions are (a,f}) = (3/5, -2/5) and 
(7,<5) = (5/3,2/3). Specifically we first fit a fine to the Ig(yA) - 
Ig(MA) distribution at each redshift, and then determine the redshift 
dependence by determining a linear fit in lg(£'(z)) to the evolution 
of the normalization constant B(z). We find no significant indica- 
tion for a redshift evolution of the slope a or 7. 
The best fit parameters and the logarithmic scatter at fixed mass. 



o-y 



ZtiiigiY,/Y(M ,))) 

N-2 



2n1/2 



(8) 



where the sum runs over all Y measurements (three projections of 
each cluster at each redshift), are given in Table[T]and Table|2] 
The two scaling relations contain the same information. While 



^ Access to the cluster simulations of sample A, including web services 
allowing to interactively produce various kinds of maps, are publicly avail- 
able via the web portal at |http://www.mpa-garching. mpg.de/HydroSims | 
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the M(Y) scaling relation is the relation of more interest for 
cosmology and is the relation used in the rest of our analysis, 
the Y(M) relation is easier to interpret if one is more used to 
thinking about clusters properties at fixed mass rather than at fixed 
Y, and we will focus the discussion of the fit results on this relation. 

The slope y of the best fit relation in samples A and B is below 
the self similar value while other simulations including cooling and 
star formation find slopes comparable to or steeper than the self 
simil ar predictions l lNagaill2006l : lBattaglia et al.ll2010l : ISehgal et all 
I2OIOI) . We find a slope in agreement with previous results if 
we only consider massive clusters with M200 > 2 x 10''* Mq//; 
("Sample B*") wh ich is identical to the mass threshold used in 
ISehgal et alj j2010l) . Projection effects may a ccou nt for some of 
the dif ference with the results of lNagail |200g) and lBattaglia et alJ 
( 120101) : these authors use spherically averaged Y measurements 
and do not include projection effects, which effectively boost the 
integrated Y signal of lower mass clusterfij and hence lower the 
slope of the scaling relation. 

After accounting for differences in the baryon fractions of different 
simulations, the normalization B of the best fit scaling relation 
for sample B* is consistent with those obtained from other 
hydrodynam ical simulations with si milar physics ( t he csf run in 
iNagail J2OO6I) and the radiative run in lBattaglia et al.l ( l2010h ). 



logarithmic scatter (0-^200, ctm^oo' "'"2500^ ~ (0.046,0.108,0.326). 

Figure[3]shows the best fit ^200-^200 scaling relation for sam- 
ple A and the distribution of the z = 1 and z = clusters, which we 
plot in the form of the SZ signal scaled for redshift evolution 

Y2oo{z) = Y2oo(z)E'""(z). (9) 

The right panel shows the distribution of the scatter around the scal- 
ing relation. 



6lgM = \g{M(Y)/M) , 



(10) 



for the full sample and subsamples. This scatter definition gives 
the logarithmic error in the mass inferred from Y measurements, 
positive scatter corresponds to clusters with Y larger than expected 
for their actual mass. At all redshifts the distribution deviates from 
lognormality with a tail at large S Ig M, causing the distribution to 
have positive skewness and kurtosis. 

The left panel of Fig.|4]shows the M200 and 7200 data from sample 
B and the best fit scaling relation. We checked by visual inspection 
that the most extreme outliers, which are all in the direction of Y 
higher than expected for the cluster mass, are indeed projection 
effects. These clusters have multiple peaks or appear otherwise 
distorted in only one or two of the three orthogonal projections, 
indicating that these are not merging systems (yet). 



The slope and normalization of the scaling relation for a 
subsample of massive clusters at z = from sample A, denoted 
as A*, are comparable to those found for the sample B*. A direct 
comparison of these numbers is complicated by the fact that slope 
and scatter of the scaling relations are mass dependent, and that 
the mass distribution within sample A does not follow the cluster 
mass function. Also sample A* consists of only 1 1 clusters, five of 
these are the most massive objects in their respective re-simulation 
region, and it is hard to assess at a precision cosmology level 
whether the non-representative environment of clusters in sample 
A affects the normalization of their scaling relation. 
The redshift evolution of the scaling relation for sample A deviates 
significantly from self similar expectations. This deviation may 
be caused by mergers: As we will discuss in detail in Sect. |4] the 
Y signal of recently merged clusters is suppressed on timescales 
of order a few Myr. As the merger rate per halo per unit time 
increases with redshift, the increasing fraction of recently merged 
clusters reduces the normalization of the scaling relation, causing 
6 to deviate from the self similar value. 



The intrinsic scatter in the spherically integrated Y parameter 
of large cluster samples has been fou nd to be close to log-normal 
JStanek et alj201(]| : lFabian et alj201ll) . However, projection effects 
due to correlated structures and diffuse large-scale structure have 
been identified as an non-negligable source of scatter and bias in the 
mass scaling relation. The non-lognormal, positively skewed dis- 
tribution of scatter in projected Compton Y parame ter in our clus- 
ter sam ple is in good agreem ent with the results of iHallman et al.l 
( I2OO7I) and lYang et alj 1201(J), who analyzed light cone/ cylindrical 
projections of the SZ effect, respectively. Based on an Edeeworth 
expansion of the mass - observable distribution, Shaw et al.l j2010l) 
find that the higher order moments do not significantly impact the 
observed cluster mass function if the product of the scatter in the 
scaling relation, ctm, and the slope of the mass function at the lim- 
iting mass of a survey is less than unity. Due to low scatter of the 
SZ scaling relation this criterion is met by all upcoming SZ exper- 
iments, suggesting that projectio n effects will be insignificant for 
cosm ological constrains (but see lShaw et al.ll2008l : lErickson et alj 
I2OIII for additional mitigation strategies). 



In the following we will focus on scaling relations within 
R200 as the M200 - Y200 relation for sample A has less scatter than 
that within Rsqq. The accretion histories at if 500 are more erratic 
than at R200 which complicates the identification of merging events 
and the interpretation of trajectories in the M-Y plane. At the 
time resolution of the simulation snapshots infalling substructures 
sometimes cross in and out of R^oo before coalescence, causing a 
series of mass jumps and mass losses in Mjoo- While it is not clear 
what the best mass definition is for a merging cluster, the scatter 
in the Mvi, - Ma relation illustrates that masses within larger radii 
are less volatile: fitting Ma as a power law in Myi, and E(z) we find 



Projection effects introduce an additive signal Fp > which scales as 
YpA cc i?^ oc Af-'^, and thus the fractional error induced by projection 
effects decreases with cluster mass 



3.2 Influence of halo concentration 

The scatter in halo concentration at fixed cluster mass has been 
identified as an important source of s catter in X-ray temperature 
lAmeglio et d]|2009i) and SZ signal (Shaw et al] 



I Yang et al.ll2009l: 
2008 : lYang et'd] 



120101) of simulated clusters. Understanding the 



I2OO 

role of halo concentration on these observables is especially im- 
portant for understanding selection biases and for the comparison 
to lensing derived cluster masses. 

The right panel of Fig. |4] shows the correlation between scatter 
in halo concentration at fixed mass and scatter in Ig 7200 at fixed 
mass for all clusters in sample B. We use the halo concentration 
measurements from Ameglio et al. (2009) derived from fitting 
NFW-profiles to the integrated mass profile over the range 
0.05 < r/ifvii < 1> and model concentration c(M2oo) with a 
power law in mass. The scatter is positively correlated with more 
concentrated clusters having higher SZ signals at fixed mass, with a 
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Figure 3. Left: Relation between mass M200 and integrated Compton 7200 parameter for tlie z = (star's) and z = 1 clusters (triangles) in sample A. The 
Compton y pai'ameter has been scaled to absorb the redshift evolution of the scaling relation in order to show the power law relation M oc ^"(1) = (}'£''''" (z))". 
The solid and dotted lines show the best fit scaling relation for sample A and its lo" error. For reference, the dashed line indicates the best fit scaling relation 
for sample B. 

Right: Distribution of residuals of the best fit scaling relation for the full sample (filled histogram) and the redshift subsamples (black/red line), and the best-fit 
Gaussian to the full distribution. The vertical dashed lines illustrate the 10% and 90% quantile for the full sample, illustrating the non lognormality of the 
scatter distribution. 
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Figure 4. Left: Relation between mass M200 and integrated Compton 7200 parameter for cluster sample B. Massive clusters with M200 > 2 X 10 Mo/ft are 

shown with filled symbols. The solid and dotted lines show the best fit scaling relation for sample B and its lo" error. For reference, the dashed line indicates 

the best fit scaling relation for sample A. The strong outliers with boosted Y signal in the low-mass range are visually identified to be caused by projection 

effects. 

Right: Residua ls of the Y-M relation at fixed mass vs. scatter in the mass - halo concentration relation at fixed mass. Concentration measurements are from 

lAmeglio et al.l I2009|^ , see text for details on the determination of c/c( A/20o)- 
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correlation coefficient of 0.30 for the full sample B and 0.68 for the 
massive subsample B*. This result is in agreement with the positive 
correlation between scatter in concentration and spectroscopic-like 
temperatu re of these clusters reported in Ameglio et al. (20091). 
Similarly, Ishaw et alj 12003) find a positive correlation between 
scatter in concentration and integrated K-parameter in halos 
from adiabatic SPH simulations and from N-body simulation in 
combination with semi-analytic gas models. On the other hand, 
lYang et al.l ([2009112010) find a negative correlation between scatter 
in concentratioiy and scatter in t emper ature and integrated SZ 
signal. As discussed in lYang et al.l ( 120 lOf) . the correlation between 
halo concentration and temperature at fixed mass depends on the 
assumed gas physics and the inclusion of radiative cooling, star 
formation and feedback ma y change the sign of the correlation. 
On the observational side, IComerford et alj ( 120101 ) find AT anti- 
correlated with Ac. However this analysis is based on a sample of 8 
strong lensing clusters and the authors note that this result vanishes 
if a different measurement for the concentration of one cluster (MS 
2137.3-2353) is used. As strong lensing selected cluster samples 
are strongly afi^ected by projection effects and are biased towards 
higher halo concentrations and X-ray luminosities than average 
clusters (e.g. Meneahetti et al. 2010 , 2011), larger. X-ray selected 
data sets like the CLASH survey dPostman et alj 1201 II) will be 
needed to observationally constrain the the correlation between 
scatter in temperature and halo concentration. 

The scatter in halo concentration at fixed mass is linked to 
the formatio n epoch of a halo wi th more concentrated halos form- 
ing earlier ('N avarro et alj Il997h . albeit with large scatter (e.g. 
iNeto et al. 2007|) which is likely due to enviromental effects (see 
also IGao & White 2007). Hence the positive correlation between 
scatter in concentration and SZ signal suggests that clusters with Y 
biased low formed more recently. 



4 SCATTER INDUCED BY MERGERS 

We now turn to a detailed analysis of the evolution of merging 
clusters around the M(Y) scaling relation fit to sample A. Figure|5] 
shows the trajectory of six massive clusters around the best fit scal- 
ing relation in the Mjoo - ^200 plane. Phases identified as mergers 
are shown in red. These examples suggest that the SZ signal lags 
behind the change in mass during extended merger events moving 
the merging clusters below the best fit s caling relation. This is simi- 
lar to the findings of lRasia et alj ( 120 llh who analyzed the evolution 
of X-ray properties of two of these clusters (g8a and gib) during 
mergers and find a time delay between mass increase and rise in 
temperature of order a few hundred mega years. We quantify the 
difference in evolution during mergers compared to the overall evo- 
lution of each cluster in the M-Y plane in Fig.|6] The open symbols 
show the logarithmic increase in mass 



and SZ signal scaled for redshift evolution 

AlgK = lghr 



AlgM = lg 



M(z = 0) 
M(z = 1) 



(11) 



(12) 
J(z=l)l 

As expected, the overall evolution from z = 1 to z = as quantified 
by the slope of the best fit linear model with zero intercept is 
consistent with the slope of the best fit scaling relation. 
The filled star symbols show the evolution of each cluster in the 
M-Y plane during merger phases only (this corresponds to the 
sum of the red line segments for each cluster in Fig. \5\ treating 
the different projections separately). The dashed red lines indicate 
the best fit slope for the relation between increase in mass and 
redshift scaled Y during mergers. This shows that the Y signal 
scaled for redshift evolution increases more slowly during mergers 
than expected from the overall scaling relation. The dashed lines 
show the best fit slope for the relation between increase in mass 
and redshift scaled Y during mergers when relaxing the merger 
criterion to include all times at which the fractional accretion rate 
is above its mean value. This illustrates that the suppression of Y 
during mergers is robust with respect to the definition of merger 
event. 

We further illustrate the connection between merging events 
and scatter in the M2oo(J'20o) scaling relation in Fig.|7] The top left 
panel shows how the clusters evolve around the scaling relation, 
giving the cumulative fraction of clusters evolving into outliers as a 
function of time, averaged over all clusters and all snapshots. Thick 
(thin) dashed-dotted or dashed lines show the fraction of clusters 
which evolve at least 10% (20%) below or above the scaling re- 
lation. For example, starting from one simulation snapshot, about 
38% of all clusters will move at least 10% below the scaling rela- 
tion within the next seven snapshots (corresponding to about one 
Gigayear), about 30% deviate at least 10% above the scaling rela- 
tion during that time period and about 35% stay within 10% scatter 
from the scaling relation. The asymmetry between these pairs of 
lines is due to the non-lognormal distribution of scatter, the thick 
lines correspond to the 24% and 80% quantile, the thin lines corre- 
spond to the 4% and 90% quantile. The top right panel shows the 
same evolution around the scaling scaling for clusters undergoing 
a merger at f = 0. Within a Gigayear after a merger, 55% of all 
clusters will go through a phase where the inferred mass is biased 
low by at least 10%, while for only 30% of these cluster the in- 
ferred mass will be biased high by more than 10% during this time. 
The bottom left panel shows the ratio of these two plots, and illus- 
trates the asymmetric evolution of mergers below the scaling rela- 
tion. The inferred mass of a recently merged cluster is about 50% 
more likely to be biased low by at least 10% and twice as likely to 
be biased low by at least 20% compared to an average cluster. 

The bottom right panel shows the cumulative fraction of 
clusters which have undergone a merger as a function of look back 
time given their current deviation from the scaling relation. This 
plot shows that 50% (75%) of all clusters with inferred masses 
biased low by at least 10% (20%) have undergone a merger within 
the last Gigayear. 



'' These authors use lg(i?20o/^50o) as a proxy for concentration, which for 
an NFW profile is a monotonically decreasing function to halo concen- 
tration. We find correlation coefficients of -0.22 (-0.47) for the scatter in 
lg(^20o/^50o) and F200 at fixed mass for sample B (B*), indicating that 
our result is robust with respect to the definition of halo concentration em- 
ployed. 



In summary our analysis shows that the SZ signal changes 
more slowly than cluster mass during mergers. This indicates that 
for a cosmological distribution of merger orbits and mass ratios the 
delay between mass accretion and heating of the ICM by shocks 
and partial virialization are more important than merger boosts. 
Hence the inferred mass of recently merged clusters tends to be 
biased low and we find that a large fraction of negative outliers are 
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Figure 5. Evolution of six massive clusters in mass and 7200. the redshift evolution scaled 7200- Offsets are added to show all clusters in one plot. We show 
three orthogonal projections for each cluster to illustrate the magnitude of projection effects. Phases identified as merging events are shown in red. The dashed 
and dotted lines show the best fit scaling relation for sample A and its lo" error 



associated with recent mergers. 

Note that throughout this section we have analyzed deviations 
from a scaling relation determined from a fit to sample A. Since the 
merger histories of this environment selected sample are not neces- 
sarily representative of a volume limited sample the calibration of 
this relation may be biased. However, the results in this section and 
the correlation between scatter in halo concentration and SZ signal 
of the volume limited sample discussed in Sect. 13.21 suggest that 
this bias would increase the normalization B and slope y at fixed 
Y . Hence such a calibration bias would downplay the asymmetric 
scatter induced by mergers that we reported in this section. This 
suggests that in a volume limited sample merging clusters may be 
less frequent, but their inferred masses could be more biased. 



5 SZ MORPHOLOGIES 

Since we found the dynamical state of clusters to be correlated 
with scatter in the M(Y) scaling relation we now test if the mor- 
phological appearance of SZ maps can be used to identify clusters 
that deviate from the scaling relation. Quantitative measures of 
the X-ray surface brightness morphology are commonly used 
to identify disturbed clu st ers, observations (e. g. iBohringer et al.l 
l20ld IOkabeetalj|2010l: lMarroneetal.ll201lh and simulations 



I Jeltema et al.l2008| - |Ventiniiiglia et al.l2008l : lB6hringer et alj20l"oh 
find the inferred masses of morphologically disturbed clusters to 
be biased low. Ventimiglia et al. ( 20 08) analyzed t h e mor phology 
of clusters from the simulation of iBorgani et al.l ( l2004f) . which 
is our sample B, and find significant correlations between the 
centroid shift, axial ratio and power ratios of the X-ray surface 
brightness distribution of these clusters and scatter in the Tx(M) 



relation. lBohringer et alj ( l2010h compared the morphology of these 
simulated clusters to observed morphologies in the REXCESS 
sample, and show that the simulated X-ray morphologies show a 
larger dynamic range and appear more disturbed during mergers. 
They trace this difference to the fact that cool cores are more 
pronounced in this simulation. 

Here we test the effectiveness of a number of morphological 
parameters, which are typically used to measure X-ray morphology 
of clusters or optical morphology of galaxies, at quantifying sub- 
structure in projected y maps. Within a circular aperture of radius 
/?2oo we compute the following quantities: 

• Asymmetry A measures substructures and differences from cir- 
cular symmetry, it is defined as the normalized difference between 
an image / and a copy R of the image rotated by 180 degree, 
^ = 2; \Ii-Ri\/ 2/ ^M where sum runs over all pixels in th e aperture, 
and th e center of the aperture is chosen to minimize A JConselicd 
l2003l) 

• Centroid shift w JMohr et alJI 19950 is another measure of the 
distribution of bright substructures based on the change of the cen- 
troid of different isophotal (iso- y) contour s. Specifically, we follow 
the implementation of Ventimi glia et al.l (|2008) and compute the 
variance of the centroid for 10 iso-y contours spaced evenly in Igy 
between the maximum and minimum of y within the aperture. 

• Concentration C We quantify the apparent concentration of 
the y distribution by the fraction of integrated Y contained within 
0.3 X R200, C = i'o.3R2oo / ^200 

• Ellipticity e = l-B/A is defined as the ratio of semi-major (A) 
and semi-minor axis (B) and is calc ulated directly from the second 
order moments of the y distribution JHashimoto et alfcOOTl) 

• Gini coefficient G measures the uniformness of pixel values 
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Figure 6. Logarithmic mass growth and increase in SZ signal scaled for cosmological evolution for all clusters in sample A. The left panel shows the evolution 
within A = 200, the right panel for A = 500. The black open symbols show the overall evolution of individual clusters between z = 1 and z = 0, the black solid 
lines are the best hnear fit with zero intercept to these points, yielding a slope of 1.62 ± 0.19 at A = 200 (1.62 ± 0.29 at A = 500), consistent with the slope of 
the best fit scaling relation. Filled, red stars show the evolution of each cluster during merger phases, the dashed lines are the best linear fit with zero intercept 
to the evolution during mergers with slope 0.94 ± 0.15 (0.95 ± 0.22). The dotted lines show the best fit slope for the evolution during mergers when the merger 
criterion is relaxed to times when the fractional accretion rate per unit redshift is larger than the mean fractional accretion rate per unit redshift. 



regardless of their spatial distribution JLotz et al.ll2004f) . It is based 
on the Lorentz curve, the rank-ordered cumulative distribution of 
pixel values. It is defined as 

^ n n 



G = 



lynin ■ 



->';l 



(13) 



= 1 ;=l 



where n is the number of pixels inside the aperture, y^ the value of 
the fth pixel, and y is the mean pixel value. The Gini coefficient of a 
uniform distribution is zero, and it is one if one pixel contains all the 
signal. It increases with the fraction of y in compact compone nts. 

• Second order brightness moment M20 JLotz et alj|2004r) : The 
total second-order moment M is the signal in each pixel j, 
weighted by the squared distance to the center of the galaxy cluster 
(xic, ^i.c), summed over all pixel inside the aperture: 



M - 



2 M, = ^ >' ((-^1.' " -'^i-^)^ + '^-^2.' " -^2.^ )^) 



(14) 



Again, the center is determined by finding (xj c, X2,c) that minimizes 
M. The second-order moment of the brightest regions measures the 
spatial distribution of bright sub clumps. M20 is defined as the nor- 
malized second order moment of the brightest 20% of the cluster's 
flux. M20 is computed from the pixels rank ordered by y. 



M,i 



log 



Hi Mi 

M 



while V yi < 0.2Y2( 



(15) 



M20 is similar to C, but it is more sensitive to the spatial distribution 



of luminous regions and is not based on any symmetry assump- 
tions. 

• Multiplicity f JLaw et al.l uOOTP is another measure of the 
amount (multiplicity) of bright substructures. Using the observed 
y distribution as a tracer of the cluster's projected mass, one can 
calculate a "potential energy" of the y distribution 



T;ictii;il — 



,=1 j=l,j^i ''J 



(16) 



where r,j is the distance between pixels / and j. This value is nor- 
malized by the most compact possible re-arrangement of the pixel 
values, i.e. a circular configuration with pixel values decreasing 
with radius. The "potential energy" of this most compact light dis- 
tribution is 






/=! j=lj*i 



yijj 



(17) 



where r'.. is the distance between pixels / and j in the most compact 

configuration. 

The multiplicity coefficient is defined as 



T = 100 X log 



"Pco 






(18) 



It is similar to A and M20, but is has a larger dynamical range than 
M20 and requires no center or symmetry assumption. 
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Figure 7. Top left: Cumulative probability for a cluster to deviate from the scaling relation by (5 Ig M as a function of time. Thick (thin) dash-dotted blue lines 
show the fraction of clusters deviating at least 0.04 (0.08) below the scaling relation, corresponding to a bias of 10% (20%) in the inferred mass. Thick (thin) 
dashed lines show the fraction of clusters deviating at least 0.04 (0.08) above the scaling relation. The black solid line show the fraction of cluster which 
deviate less that 10% from the scaling relation within a given time. In all panels error bars indicate statistical errors estimated from 100 boot strap realizations. 
Top right: The same for merging clusters. Note that extended merging events are counted as multiple mergers, effectively giving more weight to major 
mergers. 

Bottom left: Ratio of the above panels, highlighting the enhanced probability for mergers to evolve below the scaling relation compared to an average cluster. 
Bottom right: Cumulative fraction of clusters which have undergone a merger as a function of look back time and their current deviation from the scaling 
relation. 



• Power ratio P„ JBuote & Tsaill 19951) correspond to a multi- 
pole expansion of the y map inside an aperture centered on the y 
centroid. We measure the power ratio PoIPq which is related to the 
projected cluster ellipticity. 

We measure morphology at a fixed physical resolution of 17.6 
kpc/pixel and do not include any noise or observational eflFects. 

Figure [8] shows the morphology as measured by C, A, and "P 



of four massive clusters from simulation A during their evolution 
since a = 0.5. The evolution of these clusters around the M{Y) 
scaling relation is shown in Fig. |5] Vertical lines indicate the 
onset of mergers. Clusters g696a, g696c, and gib illustrate the 
expected course of a merger: As a merging object enters the 
aperture within which morphologies are computed, the clusters 
appear less symmetric (higher A), less concentrated (lower C) and 
shows more substructure (higher T). As the infalling clump sinks 
towards the cluster center and dissolves, the cluster appears less 




1.0 
0.6 
0.2 ' 

7 - 
5 



mA'^f^'F^r\- 



sA-ZV^; 



^U.. 



ts^^.f^ 



Merger induced scatter in the M-Y relation 1 1 



0696c _ 


wsSk-^C— :>^ 


v^: 




^^\^ 


\^r^-r^'^^^'^ 




0.5 0.6 0.7 0.8 0.9 1.0 



Figure 8. Evolution of morphological pai'ameters G, A, m for four massive clusters from sample A, different lines in each panel show the three orthogonal 
projections. The bottom panel shows the fractional accretion rate on a logarithmic scale, the dotted and dashed line indicate the mean accretion and the 
accretion rate threshold used to define mergers through out this analysis. Vertical lines mark the onset of mergers, i.e. the time when the fractional accretion 
rate first crosses the threshold used to define mergers. At the onset of a merger clusters appear less concentrated, more asymmetric and show more substructure. 



disturbed again. However, linking accretion history to morphology 
is comphcated by extended merger phases (g696c, gib at a > 0.8) 
with multiple infalling clumps. It is also apparent from these 
examples that fluctuation in morphology are not always linked to 
major accretion events (e.g. g8a, late time evolution of g696a). 



For a more representative distribution of dynamical states and 
morphologies, we show the distribution of scatter in the M{Y) re- 
lation and morphological parameters for all clusters in sample B in 
Fig. |9] Shaded region contain the 25% most disturbed/most elon- 
gated/least concentrated clusters. Overall, the inferred mass M{Y) 
has larger scatter for clusters with disturbed morphologies, but it 
is nearly unbiased. Splitting the cluster sample by mass shows 
that morphologically disturbed clusters with low mass (M200 < 
10'"'Mo//2, open star symbols) tend to be biased towards larger 
inferred masses, while massive clusters (M200 > 2 x IO'^Mq/Zz, 
filled red triangles) with disturbed morphologies are preferentially 
biased low in inferred mass. We quantify this trend using the Spear- 
man rank order correlation coefficient for different mass samples 
and show the correlation coefficients in Fig. |9] If the significance 
level .? of a correlation between a morphology parameter and mass 
bias is low (s > 0.01), we do not list a correlation coefficient. We 
find a significant correlation between morphology and mass bias in 
all three mass bins ( M > 2 x 10'''Mo//i,M > 10'''Mo//z,M < 
IQ^^MqIIi) for the multiplicity, concentration,M2o and asymme- 
try parameter. These different morphology parameters consistently 
show that the correlation between disturbed morphology and nega- 
tive mass bias increases with mass threshold, and the correlation co- 
efficient changes sign for the low mass clusters. For centroid shifts 
and the Gini coefficient, we only find significant correlations with 
scatter in the M{Y) relation in two mass bins, which follow the same 
pattern as just described. Power ratio P2/P0 and ellipticity are cor- 
related with mass bias only for the most massive clusters, such that 
less circular clusters tend to be biased low in mass. 
This segregation in mass, which is consistent among all morpho- 
logical parameters, suggests that a large fraction of morphologi- 
cally disturbed clusters which are biased high in inferred mass is 
caused by projection effects. The more massive clusters, which are 
less affected by projection effects, show correlations with disturbed 
morphology corresponding to a negative bias in inferred mass as 
expected from X-ray results. We expect cool cores to have a smaller 
influence on the SZ morphology than is found in X-ray, as the SZ 



signal is linear in density and less sensitive to physics in the clus- 
ter core. Projection effects due to uncorrelated large scale structure 
along the line of sight are on average more diffuse than the projec- 
tion effects from nearby structure that is included in our analysis. 
Hence we do not expect the morphology of massive clusters to be- 
come dominated by projection effects for line of sight projections 
which include all intervening structure. 

As a first step towards towards including resolution effects, we con- 
volve all projected y maps with a circular Gaussian beam with full 
width at half maximum (FWHM) of 150 kpc, and sample the maps 
at a resolution of four pixels per FWHM. For a telescope with an 1 
arc minute beam, this physical resolution is reached for a source at 
z ~ 0.15; for an experiment with beam width of about 20 arc sec- 
onds, this corresponds to z ~ 0.8. Figure [TO] shows the correlation 
between mass bias and cluster morphology as measured from these 
blurred maps for all massive clusters with M > 2 x lO'"* Mg/h from 
sample B. For this choice of beam and pixel scale, cluster morphol- 
ogy and bias in inferred mass are well correlated and resolution ef- 
fects are small. However, since this analysis is based on noise- and 
background-free y maps and a simplistic map making procedure, 
more realistic simulations are required to assess whether SZ based 
morphology can in practice be used as a proxy for the dynamical 
state of a cluster. 



6 SUMMARY AND DISCUSSION 

Using projected Compton y maps of galaxy clusters extracted from 
cosmological hydrodynamical simulations, we analyze the clusters' 
thermal SZ signal and its scaling relation with cluster mass. We 
study the detailed time evolution of a sample of 39 clusters around 
the scaling relation using simulations with outputs closely spaced 
in time. Compared to previous studies, which focused either on the 
evolution of isolated, idealized mergers or on large samples of clus- 
ters at widely spaced redshifts, this sample enables us to isolate the 
effect of merging events for a cosmologically representative distri- 
bution of merger orbits, mass ratios, and impact parameters. Our 
main results can be summarized as follows: 

(i) The best fit scaling relations to the integrated 7200 signal of 
these clusters are close to self-similar predictions and agree well 
with other simulations that include comparable gas physics. 

(ii) The scatter around these scaling relations is small (of order 
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Figure 9. Relation between scatter in the M2oo(J'20o) relation (51gM2oo and morphological parameters for all clusters from sample B measured within an 
aperture of size i?200. Open star symbols show clusters with M < IO'^Mq/A, filled circles how clusters with 10''' Mq//; < M < 2 x 10'''Mo//i, and filled 
triangles show massive clusters with M > 2 X IO'^'Mq/A. Dashed vertical lines indicated the 25% and 75% quantiles of the morphology distribution. Shaded 
regions contain the 25% of the data points which are classified as most disturbed by that morphological parameter. Numbers in the upper left or right corner 
give the Spearman rank correlation coefficient between the morphological pai'ameter and scatter in the M{Y) relation. From top to bottom these numbers are for 
mass samples M > 2 X 10 Mq/A.M > 10 Mq/Zj^M < W Mq/H. If a coiTelation is not significant (significance level > 0.01), we do not list the correlation 
coefficient. 



10% scatter in mass at fixed Y200) and it is overall well correlated 
with the scatter in halo concentration, such that more concentrated 
halos have larger Y signal at fixed mass. 

(iii) The scatter in the scaling relation deviates from a log nor- 
mal distribution and is skewed towards clusters with Y signals 
larger than expected from their mass. We find projection effects due 
to nearby structures to be an important source of this upward scat- 
ter. However, due to the small magnitude of the scatter in the mass 
scaling projection effects are not expected to be a significant con- 
tamination for cosmological constraints from SZ cluster surveys. 

(iv) Merging clusters fall below the scaling relation, such that 
their inferred masses are biased low. More quantitatively, we find 
that within a Gyr following a merger, clusters are twice as likely as 
the average cluster to undergo a phase during which their inferred 
mass is biased low by more than 10%. 

(v) We identify merging events to be a major source of down- 
ward scatter in the scaling relation: a large fraction of clusters 
whose inferred masses are biased low recently underwent a merger 
(c.f. Fig.0. 

(vi) For massive clusters, we find the morphology of SZ maps to 
be well correlated with deviations from the scaling relation. While 
the robustness of this result with respect to noise and imaging ar- 
tifacts requires further analysis, it suggests that SZ morphology 
may be useful to reduce the scatter of mass estimates, and to in- 



fer merger rates of massive halos and hence test theories of halo 
formation. 

Our analysis of the time evolution of merging events is in agree- 
ment with the conclusions drawn from earlier studies comparing 
morphologically disturbed and undisturbed clusters in cosmo- 
logica l simulations a t fixed redshif t s (e.g . Mathiesen & Evrara 



200 ll: iKravtsov et al.l I2OO6I : iNagail 120061 : Ijeltema et al.l l2008l : 



Ventimiglia et al.l l200g ). Specifically, it supports the hypothesis 
that for a cosmological distribution of merger parameters partial 
virialization and non-thermal pressure support due to mergers 
are more important than merger boosts found in simulations of 
direct collisions between mergers. For simulated clusters the 
intrinsic scatter in the scaling relation and the mass segregation 
between morphologically relaxed and disturbed clusters are 
significantly smaller than recent observational results based on 
SZ mea surements. X-ray mo rphology and weak lensing inferred 
masses JMarrone et al.l uOllh . However, as these authors note, 
the observed scatter is in agreement with the scatter expected 
in weak lensing mass measurements (Becker & Kravtsov 2010). 
Similarly, the mass segregation is enhanced by the sensitivity 
of weak lensing mass estimates to cluster triaxiality, and these 
observational constraints on the intrinsic scatter and bias in SZ 
mass estimates are limited by the accuracy of weak lensing mass 
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Figure 10. Relation between scatter in the M2oo(J20o) relation (51gM2oo and morphological pai'ameters for clusters with M > 2 X lO'^Mo/Zr from sample B, 
measured from SZ maps smoothed with a Gaussian beam with a FWHM of 150 kpc and sampled at a pixel scale of 37.5 kpc. Dashed vertical lines indicated 
the 25% and 75% quantiles of the morphology distribution. Shaded regions contain the 25% of the data points which are classified as most disturbed by that 
morphological parameter Numbers in the upper left or right comer give the Spearman rank correlation coefficient between the morphological parameter and 
scatter in the M{Y) relation. Dashed-dotted lines show the best fit linear relation. 



reconstruction. 



Further complications arise when inferring cluster masses 
from SZ observations as most Y measurements are derived from 
fitting parametric profiles (e.g. iNagai et al.l 12007 al ; lArnaud et al.l 
201(1) to the dat a which assume rad i al symmetry (but see 



zUlU) to tne dat a wmcn assuine rad i al symmetry (but see 
Plagge et alj I2OI0I ; iMarrone et all l201lh ISavers et alj I2OIII for 



alternate methods and discussions). The distorted geometry of 
merging clusters may introduce additional scatter to mass esti- 
mates derived from profile fits, but an experiment specific analysis 
of such effects is beyond the scope of this work. 

An additional limitation of our analysis is the range of 
non-gravitational physics included in the simulations. While recent 
studies show the impact of AGN-feedback on overall cluster 
profiles and scaling relations ( Sijacki et al. 2007|; IPuchwein et al.l 
[2OO8; Battaglia et al. 2010; Fabjan et al. 201 1), this mainly afi'ects 
the cluster center. Consequently, we do not expect AGN feedback 
to significantly alter the slow virialization of newly accreted 
material at larger radii, which we found to be the main source of 
scatter during merging events. In the cluster outskirts, electrons 
and ions are not in t herma l equilibrium. Rudd & Nagai ( 20091) 
and IWong & Saraziij ( 120091) show that detailed treatment of the 
multi-temperature structure of the intracluster medium leads to 
a significant suppression of electron temperature and SZ signal. 
Based on a sample of three simulated cluster, IRudd & Nagail 
( |2009() find this eifect to be especially pronounced in clusters 
undergoing major mergers. Under specific conditions, this effect 
may cause a bias of up to 5% in integrated Y , corresponding to 
an additional negative bias of about 3% in the inferred mass of 
merging clusters. 

Overall, we find that merger events cause a temporary nega- 
tive bias in inferred cluster mass of order 10% - 15%. Due to the 



increased fraction of recently merged objects at higher redshift, we 
conclude that this merger bias should be accounted for when mod- 
eling the redshift evolution in the scatter of scaling relations. 
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